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Abstract. The gap function of a domain Q C M" is 

5(C) :=<i2(A2_Ai), 

where d is the diameter of Q, and Ai and A2 are the first two positive Dirichlet 
eigenvalues of the Euclidean Laplacian on Q. It was recently shown by Andrews 
and Cluttcrbuck [l] that for any convex ft C M", 

where the infimum occurs for n = 1. On the other hand, the gap function 
on the moduli space of n-simplices behaves differently. Our first theorem is 
a compactness result for the gap function on the moduli space of n-simplices. 
Next, specializing to n = 2, our second main result proves the recent conjecture 
of Antunes-Freitas (SI: for any triangle T C M'^, 

with equality if and only if T is equilateral. 



1. Introduction 

Let Q, C M" be a convex domain. Let Ai and A2 be the first two eigenvalues of 
the Euclidean Laplacian on O with Dirichlet boundary condition. It is a classical 
result that < Ai < A2. The gap between Ai and the rest of the spectrum, 

A2 — Ai 

is known as the fundamental gap offl. The gap function 

where d is the diameter of f2. The gap function is a scale invariant: it is purely 
determined by the shape of the domain. Physically, if we consider heating the 
domain at some initial time and then keeping the boundary of the domain fixed 
at zero temperature, the fundamental gap determines the rate at which the overall 
heat in the domain vanishes as time tends to infinity. It is natural to ask the 
following question. 

How does the shape of a convex domain affect the rate at which it 
loses heat over a long period of time? 

The mathematical formulation of this question is: 

What is the relationship between the geometry of a convex domain 
n C M" and ({n)? 



Key words and phrases, spectral geometry, spectral gap, polygonal domain, Euclidean triangle, 
Dirichlet eigenvalues, Laplacian, fundamental gap conjecture. 
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M. van den Berg |17 observed that for many convex domains, the gap function is 
bounded below by a constant. For example, consider a rectangular domain R C M? , 

R = {{x,y) GM^ I 0<a;<a, 0<?/<6}. 

Using separation of variables, it is straightforward to compute that the eigenfunc- 
tions and corresponding eigenvalues of the rectangle are 

•^j.fel^;,?/) = sm sm — — L A^^fc = + — -, j,keN. 



a J \ b J 

Making the additional assumption b < a, one computes the gap function of the 
rectangle i?. 

If we then think about the gap function on all possible rectangles R, we see that 
the square uniquely maximizes the gap function with 

^(Square) = Gtt^. 

On the other hand, if a rectangle collapses to a segment, by letting b i 0, then 
f J, Stt^. An even more elementary example is the segment. The gap function on 
any (finite) segment [a, b] with a < 6 is 

Perhaps based on this intuition, Yau formulated the fundamental gap conjecture, 
in 18 which was recently proven by Andrews and Clutterbuck [l]. 



Theorem 1 (Andrews-Clutterbuck). For any convex domain in M", the gap func- 
tion is bounded below by Stt^. 

This result shows that among all convex domains, the gap function is minimized 
in dimension 1. If the gap function is restricted to a certain moduli space of convex 
domains, what are its properties? 

In this work, we focus on the gap function restricted to the moduli space of n- 
simplices and in particular, the moduli space of Euclidean triangles. Recall that an 
n-simplex AT is a set of n + 1 vectors {vq, ■ ■ ■ , w„ } in M" such that vi~vo, ■ ■ ■ ,Vn — VQ 
are linearly independent. The convex domain 
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^ tj = 1, tj > for < j < n 



defined by X is bounded with piecewise smooth boundary. For the sake of simplic- 
ity, we don't distinguish the simplex X with the domain it defines. The moduli space 
of n-simplices is the set of all similarity classes of n-simplices; it is parametrized by 
the set of n-simplices with diameter equal to one. We note that in case n — 2, this 
theorem is straightforward to deduce from the main result of 

Theorem 2. Let Y be an n~ 1 simplex for some n >2. Let {Xj}j^f^ be a sequence 
of n simplices each of which is a graph over Y . Assume the height of Xj over Y 
vanishes as j — > oo. Then ^{Xj) —> oo as j —> oo. More precisely, there is a 
constant C > depending only on n and Y such that ^{Xj) > Ch{Xj)-^/^, where 
h{Xj) is the height of Xj. 
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Since any triangle with unit diameter is a graph over the unit interval, this 
theorem implies that there exists at least one triangle which minimizes the gap 
function on the moduli space of triangles. The moduli space of triangles is the set 
of all similarity classes of triangles, which we identify with 

Wl^ e : y > 0, ^<x<l, + < 1 

where the vertices of a triangle in each similarity class are (0,0), (1,0) and {x,y). 
The following result shows that the gap function on triangular domains is more 
than twice as large as the gap function on a generic convex domain; the theorem 
was conjectured in [2j. 

Theorem 3. For any triangle T with unit diameter, 

where equality holds iff T is equilateral. 

Let us recall the famous question posed by M. Kac [o]: 

Can one hear the shape of a drum? 
The resonant tones of a domain are in bijection with the eigenvalues of the Euclidean 
Laplacian with Dirichlet boundary condition. Therefore, with a perfect ear that is 
capable of registering all tones, one can hear the spectrum, that is, the set of all 
eigenvalues. Kac's question is then mathematically reformulated as follows. 

If two domains in have the same spectrum, do the domains also 

have the same shape? 

A negative answer to Kac's question was demonstrated by Gordon, Webb and 
Wolpert [T], who showed that there exist isospectral planar domains which are 
not isomorphic. On the other hand, Durso ^ proved that if the two domains are 
triangles in M^, and they have the same spectrum, then they must be the same 
triangle. The proof uses the entire spectrum, so we can reformulate her result as 
follows. 

With a perfect ear, one can hear the shape of a triangle [s]. 

In practice, however, one does not have a perfect ear. That is, one may only detect 
a finite portion of the spectrum. Our Theorem [3] implies: 

Amongst all triangles, one can hear the shape of the equilateral 
triangle. 

In particular, the equilateral triangle can be heard with a realistic ear, because 
Theorem |3] demonstrates that the gap function uniquely distinguishes the equilat- 
eral triangle within the moduli space of all triangles. In fact, we expect that it is 
possible to distinguish triangles based on a finite number of eigenvalues. This is 
supported by numerical data in [sj , which shows that one expects that triangles are 
uniquely determined by their first three eigenvalues. 

Our work is organized as follows. The compactness result for simplices is proven 
in §2; in §3 this is refined to prove that Theorem 2 holds for all sufhciently "thin" 
triangles. In §4, we prove that the equilateral triangle is a strict local minimum 
for the gap function on the moduli space of triangles, and in §5 we determine a 
lower bound for the radius of the neighborhood in the moduli space of triangles 
on which the equilateral triangle is a strict local minimum. Finally, in §6, we 
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provide an algorithm to complete the proof of Theorem 2. Concluding remarks and 
conjectures are offered in §7. 
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2. A COMPACTNESS RESULT FOR THE GAP OF SIMPLICES 

Let US first fix the notation. The Laplace operator on M" is 

z— 1 * 

The Dirichlet (respectively, Neumann) eigenvalues of the Laplace operator are the 
real numbers A for which there exists an eigenfunction 

u £ C°°(ri) such that — Am = Am and uj^j^ = 0, (respectively, 



dn 



0). 



on 



We shall use A to denote Dirichlet eigenvalues, ^ to denote Neumann eigenvalues, 
and index the Dirichlet eigenvalues by N and the Neumann eigenvalues by OUN. The 
Dirichlet and Neumanrj^ eigenvalues, respectively, satisfy the following variational 
principles [i]. 



(2.1) 



inf 



Mi 



inf 



^^0 



M 



for fc > 1 and j > 0, where (f)j and Lpi are, respectively, eigenfunctions for Aj and 
fjLi (assuming that (/jq = and = 0). The well known property of domain 
monotonicity for the Dirichlet eigenvalues is that, if a domain C fi*, then 

In |14| , we demonstrated the following weighted variational principle for so-called 
drift Laplacians. Given a weight Junction (j), the drift Laplacian with weight (j) is 



A. 



A- V0- V. 



The Dirichlet and Neumann eigenvalues of the drift Laplacian satisfy the following 
variational principles. 



Afc = inf 

¥'eci(A/) 



IM 



M 



fJ-k 



inf 



e-^ < j < fc, ifUi = 
ip(pje~'^, -1 < j < k> 



M 



^Note that the Neumann boundary condition is automatically satisfied if no boundary condition 
is imposed in the variational principle. 
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for fc > 1 and I > 0, where ipj achieves the infimum for k ~ j (and as above (/jq = 
and ip^i = 0). Finally, throughout this paper we will use the following notations: 
for a function f{t) and fixed fc > 0, 

f{t) = 0{t'') as i ^ if there exists C,6>0 such that \ f{t)\ < Ct^ for all \t\ < S; 

fit) = ofi'^) as < ^ if lim ^ = 0. 

2.1. Proof of Theorem [2j To prove Theorem [2j we show that if a sequence of 
n-simplices collapse, there exists C > such that £,{Xj) > Ch{Xj)~'^/^ , where the 
height h of the simplex (defined in the arguments below) vanishes as j — >■ oo. For 
simplicity in notation, let us drop the subscript. We may assume that the simplex 
is defined by the points 

{p^};ucM", p°-o, 

such that 

n 

= Y.P^^'^ = 0, 1 < j < n - 1, 

where {e^}2^i are the standard basis of E". In other words, p°,...,p"~^ are 
contained in the span of {e'}"^^^. The collapse is described by 

In fact, we may assume without loss of generality that the simplex is contained in 
the set of points 

I X e M"| X = ^Xfce*^, a;„>ol. 



fe=i 



Then, for any point 



the height of x 



xeX, x^^Xfce'', 



k = l 



/i(x) x„. 

The height of the simplex itself is defined to be 

h = h{X) :=Mp")=K- 

Since the simplex collapses, we assume in the remaining arguments that h < 0.1. 

Let Xi, i = 1,2, be the first and second Dirichlet eigenvalues of X with corre- 
sponding eigenfunctions ipi such that = 1. In the following claim, we demon- 
strate the quantitative estimate that at least 90 % of the mass of the eigenfunctions 
(pi and (j)2 is contained in a cylinder around p" intersected with X. We call this 
estimate "cutting corners" because it shows that we may "cut off the corners" and 
use the cylinder to estimate the gap. Let 

n-l 
i=l 

and let B^~^/^ (p) be the {n—1) dimensional ball in the space spanned by e-'^, • • • , e""-'^. 
The constant c will be chosen later. We define U to be the intersection of the cylin- 
der with base B'^~y^(j>) and height h with X, 
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where Ih is the interval of length h. Let 

V:=X\U, 

and let 



Claim: There exists a constant A which depends only on n and Y such that if 

/ 1 X 3/2 

c> A and h < I — 

then 



Proof of Claim: We shall begin by assuming 

h < 



I \ 3/2 



2c 

which guarantees 

l_c/^2/3> i. 

- 2 

By definition of the simplex as the convex hull of its defining points, since p*^, . . . , p"~^ 
are contained in the span of e^, . . . ,e"~^, the diameter of the simplex is 1, and 
h < 0.1, it follows that 

(2.2) /i(x) < h{l - ch^'^), Vx e V. 



By the one dimensional Poincare inequality and since Ju (l^i = ^ ^ Jy 4>'i 



(2.3) ^■^)?(^-x/-j n^(i:c.^/.)^ i^?- 

On the other hand, X contains a cylinder 

[0,/i(l-/i2/3)] X /i2/3y, 

where h'^^^Y is the base scaled by h^^^. One computes explicitly 

^2 



(2.4) A2(E) 



C2 



where C2 is the second Dirichlet eigenvalue of Y. Consequently, (2.3) and (2.4) 
imply that for i = 1,2, 

^2 



/l2 



~ ^ h^il-chV^yjy'^^ ^'^^^ - /.2(l-/,2/3)2 + 



which shows that 
^2 



where the final inequality follows since /i < j^. On the other hand. 



2cJ - 2 - (l-c/l2/3)2 
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SO for 1 = 1,2, 
(2.5) 

Therefore, letting 



< 



C2 + 37r^ 
2c7r2 



, 5(C2 +37r2) , . , , / 1 

A = ^ 5 ^ , then c> A and < — 

TT^ V 2c 



C2 + 
2c7r2 

3/2 



10 



□ 

Consider the so-called "drift Laplacian" /S.jj on U with respect to the weight 
ftmction / — — 21og(^i, 

A[/ := A + 2Vlog(/)iV. 
Let fjL be the first non-zero Neumann eigenvalue of A;/ on J7, and let 



Then, t/j satisfies 
Let 

Then 
(2.6) 



A^} + 2V \og(t)iVi) = -(A2 - Ai)V'. 



■0 :=?/' — a. 



0^0^ = 0. 



Thus, by the weighted variational principle, since tp satisfies (2.6), 



(2.7) 

We have 
(2.8) 



< 



|V^|^0f< / IWr^f = (A2-A1) / 0^0f = A2-A 



X 



X 



Using the claim we have. 



/,2 J,2 



X 



/,2 J,2 



>io- 



Since (j)i4>2 = 0, 









/ 0102 




/ 0102 


Ju 







so by the Cauchy inequality. 



SuM2 ^ ^ 1 



9/10 



Thus, 
(2.9) 



72 ,2 9 1 1 
-0 0? > > -• 

^ ^1 - 10 9 2 



Putting together (|2J|), (|2^, and ([2^, we have 
(2.10) M<2(A2-Ai). 
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Since fj, — \2{U) — Xi{U), and U is convex, by [l] 
We estimate that 

d{Uf < {2ch^'y + h\ 
This estimate for d{U) together with ( 2.10[ ) and (2.11) give 

37^2 



Fixing c, (2.12) demonstrates that S,{X) > Ch ^/^^ — )• cx) as /i — 0, for a constant 
C which depends only on n and the base of the simplex Y . □ 

3. Theorem 2 is true for s/iori triangles 

Refining estimates from the proof of Theorem 1, we demonstrate that if a triangle 
is sufficiently "short," its fundamental gap is strictly larger than 647r^/9. 

Proposition 1. Let T be a triangle with vertices (0,0), (1,0), and {xQ,h), where 
h < 0.005, and 0.5 < a;o < 1. Let Ai and A2 be the first two Dirichlet eigenvalues 
of T. Then, 

64^2 

Proof: Define 

U {{x, y) eT:xQ~ ch^/^ < x < xq + ch^^^}, V := T - U, 

where the constant c will be specified later. The main idea, as in the proof of 
Theorem 2, is that A2 — Ai is well approximated by the first positive Neumann 
eigenvalue of U . Assume the eigenfunctions for A; satisfy 

^02^1, i = l,2, 

and let 



/3:=max^ / } " - "T^. /-^ 



Noting that 



(/-a)02=O, 



the weighted variational principle for the first positive Neumann eigenvalue fJ.(U) 
of the drift Laplacian Ajj with weight function — 21og0i gives 

- w - luif - ^y^i luif ~ ■ 

We compute the denominator 

2j,2 / J.2 „2 / / j2 \ ^ 1 a 2 



if - aY4>i = J^q^i- a' \^J^ cpf j > 1 - P - a' 

where we have used the definition of /3 with 02 — ^ ~ Iv '^2' together with 
/[/ 01 < /t 01 = 1- So, we have 

A2 — Ai 



Mt^)<. o 2- 
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By Corollary 1 of ^ and Corollary 1.4 of [1, 



<p{uy 

which implies 

A2- Ai > {l-p-a^)- 



(p{uy 

Since (pi and 02 are £^ orthogonal, 

which by the Cauchy Schwarz inequality and definition of j3 gives 
Consequently, 

Proceeding by contradiction, we assume 

.o^N , ^ 647r2 , ^ 647r2 

(3.2) A2 - Ai < — =^ A2 < Ai + — . 

By trigonometry, T contains a rectangle 

i?= [0,/l2/3] X [Q,h-h'''\ 

By domain monotonicity, 

Ai < Ai(i?) = ^ + — . 

/l2(l-/lf)2 /il 

The height of V is at most 

hil- ch"/^) , 

since a;o S [0.5, 1]. By the one dimensional Poincare inequality for i = 1,2, 



Ju'^' ' h^[l~chiY Jv 
Since 

647r2 ^ 647r2 
Ai < A2 < Ai + — — < Ai(i?) + 



9 - ^ ^ 9 ' 
by definition of /3, since for either i = 1 or i — 2, Ju 4>f — 1 — /3 , 

-nrfl - /3) H < A2 < 5 ^ r H , 

so we have 

-c^h*/^ + 2ch^/\ 1 ,2/3 64/i2 ^ 



(1 - c/l2/3)2 - (1 - /i2/3)2 9 

Since h < 0.1, we have 

f (1 - c/l2/3)2 / 64/l2 
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which simphfies to 
(3.3) 



< 



c(2 - c/i2/3) ) y''^ 9 



Since we assume h < 0.005, then for any c < 34, ch'^^^ < 1. In particular, fixing 
c = 10, we compute that (3.3) gives for any h < 0.005, 

/? < 0.12. 



Then, 



13' 



> 0.86, 



so (3.1 ) becomes 



Since 



Aa - Ai > 



(l-/?)2 

3 * 0.86 * TT^ 



we compute that for any h < 0.005, 

647r2 



A2 - Ai > 



9 ' 



which contradicts (3.2 1. Thus ii h < 0.005, then we have ^{T) > 



647r^ 
9 



□ 



4. The equilateral triangle is a strict local gap minimizer 

The main result of this section demonstrates that the equilateral triangle is a 
strict local minimum for the gap function on the moduli space of triangles. In the 
proof, we consider all possible linear deformations of the equilateral triangle and 
demonstrate that in any direction, such a deformation strictly increases the gap 
function. 




Figure 1 . Linear deformation of a triangle. 
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4.1. Linear deformation theory. Let T be a triangle with vertices (0,0), (1,0), 
and z = (j, k), and side lengths A < B < 1. Consider a deformation to the triangle 
T{t) which has vertices (0,0), (1,0), and z + ix, where 



(a, 6), 



1, t>0. 



The direction of the deformation is given by (a, b), while the magnitude is given by 
t > 0. The linear transformation which maps the triangle T to the triangle T(t) is 
represented by the matrix 

[1 



ta 

.0 i + f. 



We may view the linear transformation T T{t) as a change of the (Euclidean) 
Riemannian metric on M^. In other words, T{t) is isomorphic to T with the metric. 



9 = {dxf 



We compute 



Thus. 



9-' = 



2—axdy 
k 



det g 



k^+2bkt+t^ 
(fc+6t)2 
■tak 




{k+tby 



tb 

+ 1 

— tak 
(k+tbY 
k^ 

(fe+t6)2. 



B 
D 



If the eigenvalues of the original triangle and the deformation triangle are respec- 
tively Ai and Xi{t), then they satisfy 



(4.1) 7_Ai < Xi{t) < 7+Ai, 

where 7± are the eigenvalues of g~^. It follows that 

|Ai(t)-Ai| < (7+-7_)Ai. 

We compute 

A + D , J{A - Dy + 4B2 

7+ = lb . 

Substituting the values of A, B, and D gives in general 

^4^2 + ^2 + 4hkt 

^+-^-- * {k + tbY ' 

The relationship between integration over T{t) and T differs by a linear factor. 



Itu) ( ^ ) It 



lT{t) 

where throughout this paper, integration is with respect to the standard Lcbesgue 
measure dxdy on M."^. The Laplace-Beltrami operator associated to a Riemannian 
metric (in dimension n) is 



1 " 
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SO one computes the Laplacian for the deformation metric g is 

Henceforth we shaU use Aq = 9^ + dy for the Euchdean Laplacian, A for the 
Laplacian associated to a deformation metric, and V/ = {fx,fy) the gradient 
(with respect to the standard Euclidean metric). In §6, we shall use this linear 
deformation theory to complete the proof of Theorem 2. Presently, we specialize 
to the equilateral triangle which we call T, and whose vertices are 

(0,0), (1,0), and (l^^y 

A triangle obtained by a linear deformation of T, with vertices 

(0,0), (1,0), and U,^\+t{a,b), 



is equivalent to T with the metric 



9 = (dxf + -^dxdy + M 1 + ) + ) {dyf. 



ita , , I ( ^ 2tbV "^t^a^^ 



We compute, 

det fi( = ( 1 + 

and 



^/3 



\i{t) -Xi , , ^ 4^3 + 2VSbt + 

(4.2) ai := — — , a, < ^ . 

^ ' t ' I " - (\/3 + 2t6)2 

The associated Laplace operator 



Let 



1 / / / 2tbY 4t2o2\ „ 4to^ ^ ^2 



1 + 



2tb 
V3j 



then 

(4.3) A = dl + d^+tL = Ao + tLi+t'^L2, 

where 

^^■^^ {V^f2tbr L,{t = 0) = ^{-adxdy-bdl). 

and 

(4.5) L2= ^ ^-'^^ '^2^2^ 



{V3 + 2tb)^^ y^ 
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By the variational principle since Ai is smooth, the first eigenvalue for the T{t), 
which we write as Xi{t) satisfies 

(4.6) Xiit) = Xi+t\i+t^Oi, Ai = - / 0i(Li|t=o)0i, 

Jt 

where (f>i is an eigenfunction for Ai with unit norm. If A2 is simple, then we also 
have 

X2it) = Ai + tX2 + f02, X2 = - [ 02 (Ll|t=o) 02, 

Jt 

where 02 is an eigenfunction for A2 with unit norm. 

In general, A2 is not differentiable because the second eigenspace may have di- 
mension 2; this is the case for the equilateral triangle. Nonetheless, we may use the 
variational principle to show that the equilateral triangle is a strict local minimum 
for the gap function restricted to the moduli space of triangles. 

Theorem 4. The equilateral triangle is a strict local minimum for the gap function 
on the moduli space of triangles. 

We will prove the theorem by applying the following proposition together with 
explicit calculations for the eigenvalues and eigenfunctions of the equilateral trian- 
gle. 

Proposition 2. For any deformation of the equilateral triangle which preserves 
diameter, for the corresponding Li, 

■2(ii|t=o)02 - / 0i(ii|t=o)0i ) >0, 

T JT J 

for all eigenfunctions (j>i for Xi{T), i = 1,2, with unit norm. 

Proof that Proposition [2] implies Theorem |4} Let 0i and 02 be eigen- 
functions for the first two Dirichlet eigenvalues Ai and A2, respectively, for the 
equilateral triangle T. Assume the eigenfunctions have unit norm on T. Let /i 
and /2 be eigenfunctions for the first two Dirichlet eigenvalues of T{t). Let 

where above and indeed throughout this paper, integration is over the equilateral 
triangle T unless otherwise indicated. For simplicity, in this section we shall use Li 
to denote Li\t^Q. Since ^ is a linear transformation from T to T(t), fi and /2 are 
orthogonal with respect to dxdy, so by the convergence of /i — > 0i, 

tA'7\ //2(/l-0l) „,,s , „ 

(4.7) e = p- — — 0[tj as t — > 0. 

J <?^'i/i 

Note that 

^ =0, 



Jif2+ef, 



so by the variational principle, 

, ^ -/(/2 + £/l)Ao(/2+£/i) 
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Since these functions are uniformly bounded in C'' for any fixed k, the Laplace 
operator on T{t) 

A = Ao+tLi+0{t^). 

Then, 



(/2 + e/i)(A - tLi)(/2 + eh) + 0(1^) 



{h + eh)Mf2+sh)^- 

= \2{t)+t J hLih + 0{t'). 

Consequently, 

A2 < X2{t)+t J hLih+Oit"). 
Since Ai is diffcrcntiablc with Ai = — J 
(4.9) Ai(t) = Ai + tAi + 0{f) =^ Ai = Ai(i) + 1 / (/-iLi^i + Oit"). 



Then, and (|4J|) imply 
(4.10) 



A2-A1 < [A2(i) - Ai(t)] +i 
Since the deformation preserves diameter, we may re- write (4.10) as 

-o{t'). 



0{t^). 



(4.11) 



i{T)<^{T{t)) + t 



We can always construct a sequence of eigenfunctions /2 which converge in to 
some eigenfunction (/)2 for A2 with J 02 = 1- Consequently, 



m<am)+t 

Since for all 02, 



'2-^102 



01-^101 



0(^2). 



^2^1<?'2 



i>iLi0i < 0, 



we have £,{T) < ^{T{t)) for all t sufficiently small. Finally, we note that we need 
only consider deformations in directions which preserve the diameter because the 
gap function is scale invariant. We have therefore reduced the theorem to verifying 
explicit calculations involving the eigenfunctions and eigenvalues of the equilateral 
triangle. 

4.2. Eigenfunctions and eigenvalues of the equilateral triangle. In 1852, 
Lame computed the eigenfunctions and eigenvalues of the equilateral triangle by 
(real) analytically extending the eigenfunctions to the plane using the symmetry 
of the equilateral triangle 12 , 10 , 11 . The eigenvalues are given by the general 
formula 

167r2 



(4.12) 



A 



27 



(to +n —mn), m,n£'Z,, 



such that m and n satisfy 

(4.13) TO + rt = mod 3, to ^ 2n, n ^ 2to and to ^ — n. 
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The (real) eigenfunctions are 



± cos 



(2n — m)y 



{2n — m)y 



and ± sin 

(m,n) 

where the sum alternates over the orbit 
(—71, rn — rt) I— > (— n, — m) i— > (n — m, — rn) i— > (rt — to, n) i— > (to, n) i— )■ (to, m — n), 



such that {m,n) satisfy (4.12 4.13) for the eigenvalue 

16^2 



A 



27 



(to + n — mn). 



4.2.1. The first eigenspace of the equilateral triangle. The first Dirichlet eigenvalue 
of the equilateral triangle with side lengths one is given by (4.12) with m = and 
n — (or with to = n = 3), 

16^ 
3 ■ 



(4.14) Ai 
Since the first Dirichlet eigenvalue is always simple, it follows from Corollary 2 of 



15 



that the first normalized eigenfunction of the equilateral triangle T is 



01 (a;, y) 



33/4 

Proposition 3 



2V2 ( . (A^y 
I sin — 1= 



sin \ 2ti \ X 



sin 27r a; 



Viy 



(t>iix,y) 



2V2 . / 27r\/3y\ . 



33/4 



sm \ TT \ X 



/3y 



sm \ TT \ X 



Proof The standard angle addition and subtraction identities for sine and cosine 
show that, 



sm 



/47iV3j/ 



2cos27ra;sin 2t: 



V3y 



= 2 sin 27r 



Viy^ 



27^^/3y ^ 
cos — cos Zttx 



'27rV3y' 



sm \ n \ X 



sm \ n \ X 



V3y^ 



The last inequality follows from the identity 

1 



cos a — cos /3 = — 2 sin - (a + /3) sin - (a — /3) 



□ 



We compute 



327r4 



8^ 
3 



vv 



l)x{4>l)y = 0. 

327r4 



33/4 
y27r' 
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l)xx{(f>l)xy = ^ {(t)l)yy{(f>l)x' 



4.2.2. The second eigenspace of the equilateral triangle. The second Dirichlet eigen- 



value of the equilateral triangle is given by (4.12) with m = 1 and n = 5 (or with 
m = — 1 and n = 4), 



A, 



11271^ 

9 



This eigenspace has dimension two. An orthonormal basis of eigenf unctions is 
given by 



and 

v{x,y) 



2tt 



{5x - VSy) 



27r 



(5a; + ^/3^/) + cos ^{-x + 3V5y) 



31 V -cos^(-a;-3\/32/) + cos^(-4a;-2y32/)-cos^(-4a; + 2\/3y) 



2 

3t 



2tt I 



r (5x - V3y) - sin ^ (5a; + ^3?;) 

sin ^(-a; - 3V3y) + sin ^(-4.t - 2^/3?/) - sin ^ (-4a; + 2V3y) 



- sm ^[—x + 3\/3y) 



The following calculations will play a key role in the proof of Theorem 2. 



656 1 567r^ 



V,, = -- 



800 9 
6561x/3 



-6561^3 



800 



UxVx 



y^Vy -— j U^yUy 

6561^3 



800 



6561 



800 ' ' ' ^ 
51037r2 156871^ 



800 
6561 567r2 



^« 200 81 ' 

77r2(-59049 + 448007r2) 
' 5400 

51037r2 15687r4 



800 9 
51037r2 15687r4 



yy 



200 81 

5103Vl37r2 
~ 200 ' 



40 27 

77r2(59049 -|- 448007r2) 
5400 

51037r2 156871^ 



XX ^xx 



40 27 

15309\/37r2 
200 ■ 



_ 5103%/37r2 
"yy^-yy ~ ^ • 



4.2.3. The third and higher eigenspaces of the equilateral triangle. The third (dis- 
tinct) eigenvalue is given by (4.12) with m ^ n ^ 6, 



A., 



647r2 



The eigenspace has dimension one. The eigenfunction is 

2 / 27r 27r 27r A 

A{x,y) = — 2sin— (6a; + 4%/32/) - 2 sin — (6x -I- 2%/3y) - 2sin — (2%/37rj/) . 

34 \ 3 3 3 / 
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4.3. Proof of Proposition [2[ Any real eigenfunction if of A2 is a linear combi- 
nation, 

if = au + (3v, + = 1, 

and we may assume without loss of generality that a > 0. We compute 

567r2 [9 .9 , /;t .1 6561 



Q,2 _ ^2 _^ 2\/3a/3 



800 



and 

6561\/3 



o2 2 2a/3 



800 

As previously observed, we need only consider those deformations in directions 
which preserve diameter, and by symmetry, we need only consider those directions 
9 with cos6' > 0. These are deformations in directions 6 g [— 7r/2, — 7r/6], so the 
direction vector (a, b) satisfies + = 1 , a > 0, and a + ^/Sb < 0. We compute 
the minimum of 

f2bXi ^b f ^ 2 4a 



/^256007r2 + (a2-/32 + 2V3a/3)59049\ 6561 /^0 2 2a/3\ 

'0 — — p — Q H ^ a. 



y 1800^/3 y 200 V \/3 



Since (a, 6) is a unit vector, a > 0, and 6 < — it follows that b = —\fY^^~a} . So, 
we determine the minimum of 

(256007r2 + {o? - 13^ + 2\/3a/3)59049) + 59049V3(a2 - ~ 

1800^3 

subject to the constraints 

/3 

Introducing the polar coordinates, 

cos(t) := a, sin(t) := /3, 
we compute that / is minimized for 

t — a — , 

2' 2 ' 

and the minimum is 

256007r2 - 236196 „ 
> 2.64 > 0. 



3600^3 



□ 



Remark 1. In the proof of the theorem, we have shown that for any triangle T e 9Jl 
with vertices (0,0), (1,0), and {x,y) where (x - 1/2)^ + {y - VS/2f < t^, 

A2(r) - Ai(r) > ^ + t(2.64) + 0{t^). 
In the arguments below, we use the calculations in §2 to precisely estimate the error 
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5. Theorem 2 is true for almost equilateral triangles 

The following proposition is the last step we need to reduce the proof of Theorem 
2 to finitely many numerical calculations. 

Proposition 4. Let T he any triangle with vertices (0,0), (1,0) and {x,y) such 
that 

x^+y^<l, {x-lf + y^<l, 



Then, the gap function ^ (T) satisfies 



< 0.0004. 



647r2 



Proof: We shall again use Ai, i = 1, 2 for the eigenvalues of the equilateral 
triangle and Xi{t), i — 1, 2, for the corresponding eigenvalues of a triangle T{t) 
which satisfies the hypotheses of the proposition. 

Since Ai is differentiable, we have 

Xi{t) = Xi+tXi + Oi{t^). 

By the variational principle, 



Xiit) < 



Ith) '^1 



Xi{t) < - / 0iA0i, 



which simplifies to 



since integration over T and T(t) differ by linear factors which cancel in the nu- 
merator and denominator, and J (l)^ — 1. We then have 



Xiit) < - / (/.i(Ao + <L)0i, 



where L = Li + L2 is defined in (4.4 4.5 1. So, wc compute directly 



Xiit) <Xi-tJ cbi (ii|t=o) (/.I - i / (/)! (L - Li|t=o) ^1- 
Thus, we have made explicit 



\Oi{t')\ < 



We have 

L — Li \t=o 



y 

(V3 + 2i&)2 



t / (^1 (L-Li|t^o)0i 
-4x/3 



(\/3 + 2i&)2 73 



+ {ad^dy + bdl) 



We estimate using the calculations for the first eigenfunction of the equilateral 
triangle and + 6^ = 1^ 



\Oi{t')\ <t 



4^3 



(V3 + 2i6)2 



\b\ / (0i),J+t 



2 Ma' - b' ) 

{V3 + 2tby 
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( 16^3 ■ 



m Stt^ 4 87r2 



\V^{V^ + 2tbY 3 (V3 + 2t6)2 3 y 
which for t < 0.0004 gives 

(5.1) \Oi{t^)\<t'^ (175.95). 

5.1. Estimates for the second eigenspace. Since A2 of the equilateral triangle 
is not differentiable, the estimates for its error term require a bit more work. The 
main idea is to expand the first two eigenfunctions for the linearly-deformed triangle 
using the orthonormal basis of eigenfunctions for the equilateral triangle. We then 
use the Poincare inequality and our explicit calculations for the eigenfunctions of 
the equilateral triangle to estimate the error. 

5.1.1. The first eigenfunction of the linearly deformed triangle. Our eventual goal 
is to estimate X2{t) from below. To accomplish this, we require not only estimates 
for the second eigenspace of the linearly deformed triangle, T{t), but also estimates 
for its first eigenfunction. Let / be the first eigenfunction of T{t) and write 

/ = tg, ai := — ^ for i = 1, 2, 

with 

(5.2) / = ^' / '^^^ = 

As usual, integration is over the equilateral triangle T with respect to the standard 
measure dxdy, and we use 1 1 • 1 1 to denote the norm over T . Since we assume 
t < 0.0004, by (|42]) 

(5.3) \X,{t) - X,\ ^ t\a,\ <t(2.32)A,. 
We compute 

(5.4) I fLf= I ^iLh +2t [ gL(l>, + f gLg. 



Since / ipig — 0, by the variational principle (2.1 ) 

which gives the Poincare inequality for 

(5.5) ||g||<^||V5||. 

V ^2 

To estimate ||V(?||, we use the definition of / and g to compute 

(5.6) (Ao + Ai)5 = -ai/ - L01 - iLg. 
By definition of g and ( |5.2[ ) , 

(5.7) I f9= I {cl>i+tg)g^t [ g\ 



We compute using integration by parts and then substituting (5.6 5.7 1 

J g{Ao + Xi)g^- j | Vgp + Ai j g'' ^ ~tai J 9^ - J gLcbi - t J gLg, 
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which gives 
(5.8) 



By definition of L, and since + 6^ = 1, for any function i/) which vanishes on dT, 
we have 



l|Vg|P - (Ai + ait) / .9' + / gL^i + t / gLg 



tpLip 



4^3 ,, , ,,,, , ,, 4 &K/3 + 4t62 
- (73-2*161)2'"^"""^^" (y3_2i|6|)2' 



and since we always have \\ipx\\ < 11'^ "^IL IIV'j/ll ^ 11'^ V"!!; ^-iid — 1 < 6 < 1, we have 
8V^ + 4i 



(5.9) 



< 



Applying this to we have 

gLg 



||Vi/)|r for any ■0 which vanishes on dT. 



< i^||V.|P. 



(y3-2i)2 



By the Poincare inequahty for g ( |5.5[ ), ( |5.8[ ) with the Cauchy inequahty, and the 
above estimate, we have 



||V,||2. 11^511 



I|i01 



A2 ' ^A2 " ' (V3-2t)2' 

which gives 

' 8Vl5 + 4t Y ~' 
(\/3-2i)2 ) 

We assume i < 0.0004, so substituting the estimate (|5.3[) for ai, we have 



I|V5||< 



A2 - Ai - \ai\t 
A2 



I|V5||< 



A2 - Ai - O.OO232A1 
A^ 



0.0046309 



1 1^01 1 



/A, 



Using the values for Ai and A2, we have 

(5.10) ||Vg|| <O.15948||L0i||. 

Next we estimate For these calculations it is convenient to drop the 

subscript and write simply 0. We have 

Ata^ , 4\/3a 



L(j) = 



(\/3 + 2i6)2 



4\/36 + 4t62 
(7f+2t6)2'^"^ ~ (V3 + 2t6)2'^^^' 



Recalling 



and a2 + 52 — 2, we compute 
,n2 16i^ M , , 



48 



0, 



48 + 32V3t + 16^2 , 



(V3-2t|6|)4"'^""" ^(V3-2t|6|)4"'^"^" - (73-2^161)4 
Assuming t < 0.0004 and substituting the value of the integrals, we have 



Ili^ilP < (1.7861)t- 



3277" 



-5.3581- 



327r4 



-5.3581- 



327r4 



-6.1870t- 



327r'' 



-1.7861t^ 



,3277" 
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Assuming t < 0.0004, we have the estimate for 

(5.11) ||L<^i|| < 86.194 
and this gives us the approximation for HV^H, 

IIV^II < 0.159481 |i^i!>| I < 13.747. 
Moreover, we have the estimate for H^H, 

(5.12) IWI<^<e«<1.2404. 

V-^2 \/ iizn 

5.1.2. The second eigenspace. Let F be an cigcnfunction in the second cigcnspace of 
T{t), and assume | |F| | = 1, where as usual, the £^ norm is taken over the equilateral 
triangle T. Expanding F in terms of the eigenfunctions of the equilateral triangle, 

F = <p + A<i)i+ tG, 

where (fi is an eigenfunction for A2, and G satisfies 

J G^i =0, J G4>2 = V eigenfunction ^2 for A2. 

Then, wc have 

(5.13) j GF = tj G"" and j Fct>i=A = j F{f - tg) = -t J Fg, 
so by the Cauchy inequality 

A^<t^\\g\\^\\F\\\ 

By definition of F and G, 

J F^ = l + A'^ + t'^\\G\f. 

Combining these, we have 

A^ < = i'||5ll'(l + + t^\\G\f), 

so we obtain the estimate for A, 

.5 ^2 < miHl + t'WGW^) ^ < ^||g||(l+t||G||) 

Since G is orthogonal to the first two eigenspaces, the variational principle for A3 
gives 

which implies the Poincare inequality for G, 
(5.15) l|G||^<M^. 

We estimate G in the same spirit as g. We compute 

-A(l3i + F 



(Ao + A2)G = (Ao + A2) 



since 

(Ao + A2)<^ = 0. 



t 
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So, 



(Ao + A2)G = - (AiA0i - + (A - tL + A2(<) - ta2)i^) . 



(Ai - A2)A</>i 



(L + a2)F. 



To estimate HVGH and hence ||G|| by the Poincare mequahty (5.15), we use the 
above calculation together with integration by parts (as we did with g), 

which since G is orthogonal to becomes 

GLF - a2 I GF. 



t 



-- j G{L + a2)F 



By definition of F and (5.131, this is 

- j GL{A(j)i+Lp + tG)-ta2 J G^ = -aJ GL(t)i- J GLtp-t J GLG-ta2 J G^ 
On the other hand, integrating by parts gives 

J G{Ao + X2)G = - J IVGP + A2 J G^ 
Combining this with the above calculation, we have 

- J |VGp + X2 J G^ = -A J GL(f>i - J GLif - t J GLG - ta2 J G^. 
The estimate for L (5.9) and the Cauchy inequality imply 



||VG|p<A2||G|p + ||G||(|A|||L0i|| + ||L^||)+t 



8%/3 + 4t 
(V3 - 2i)2 



\\WG\\' + t\a2\\\G\\\ 



By the Poincare inequality for G (|5.15|) 
A2 



||VG||^<-^||VG||^ 
A3 

This gives the estimate 



|VG| 



i\A\\\LM + \\L^\\)+t 



8V3 + 4t 
(V3 - 2t)2 



||VG||^+^^||VG|| 



||VG||<(l-^^±^-t^«^ + 4^ 



A. 



(V3 - 2<)2 



nA\\\Lct>^\\ + \\L^\\ 



/A, 



and 



l|G||< 



1 



/a; 



^ _ A2 + t|Q:2| _ ^ 



8\/3 + 4t 



A3 ^\{V3-2t)y 
Expanding and simplifying we have 

{V3-2t)\\A\\\Lq^^\\ + \\L^\\) 



f\A\\\L^^\\ + \\I^ 



\\G\\< 



(V3 - 2t)2(A3 - A2 - t\a2\) - tX^iSVi + 4t) ' 
Recalling the estimate ( |5.14[ ) for A, 

(\/3-2i)2 



l|G|l< 



{V3 ~ 2i)2(A3 - A2 - t\a2\) ~ t\^{%^/i + 4t) 
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Collecting the ||G|| terms, 



\G\\ 1- 



(^/3 - 2t)2(A3 - A2 - t\a2\) - t\s{SV^ + At) .Jl~fi\\^ 



< (^/3-2^)2 / ||L.^i||^||g|| A 

- (x/3-2i)2(A3-A2-t|a2|)-a3(8x/3 + 4i) \^l-t^\\g\\^ ^"J ' 

This gives the estimate from above for 
(5.16) 



(^/3 - 2t)2 it\\L4>M9\\ + ^l-t^\\g\V\\L^\\ 



((V3 - 2i)2(A3 - A2 - t\a2\) - tA3(8\/3 + 4f)) v/l-^'llfflP -t^{Vi~ 2t)^\\L^M9\\ 

At this point, we may substitute estimates for every term except Hif^H, which we 
now estimate. By definition of L, 

Ata? A^/Tia A^/2,h + Ab"^ 

~ (^ + 2*6)2"^"" ^ (7f+2i6j2 '^''^ ~ (7f+2t6j2 
By the triangle inequality for C? and since t < 0.001, 
(5.17) \\L^\\ < 0.0013365||^,,|| +2.3148||^^j,|| +3.6512||^yy||. 

Since tp is an jC^ orthonormal cigenfunction for A2, 

ip = au + Pv, o? ^ fP' = 1. 
By our calculations for the second eigenspace of the equilateral triangle, 

/ „77r2(-59049 + 448007r2) _ 77r2 (54049 + 448007r2) „ „ 15309^37:2 
V 5400 5400 ^ 200 

/ „ / 51037r2 I5687HA „^ /51037r2 15687r4\ "siOSv^ 

/ „ /5 1037r2 15687r4\ / 51037r2 15687r4\ ^^SlOSv^ 

ll^^ll = [-A^ + ^ j + + ^ j + 2"^— 1(^- 

We introduce the polar coordinates 

cos(t) := a, sin(f) := (5. 
The expressions simplify a bit, 



/l5687r4 J53 097r2 . J5309^7r2 

M = y - cos(2t)^^ - sm(2t)^^. 

,, /l5687r4 , 5103772 ""siOSy^ 

<^x« V cos(2i) — — sin(2t) . 

iiv-xj/iiy ^ 200 ^ ' 200 

/l5687r4 ~~5 1037r2 "TTsiOSv^ 
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Using these calculations and (5.171, we estimate 1 1 1 1 by determining the maximum 
of 

^^^.oo^. /l5687r4 ~~~15309^^2 ~~~1530V^ 

0.0013365\/ cos 2t) sin(2t 

V 27 ^ ^ 200 ^ ' 200 



oo .o /l5687r4 ,„, 51037r2 . , SIOSn/Stt^ 

+2.3148^— - cos(20— - 3m(2t)^^ 

+3.6512^ + cos(2t)^^ + sin(20 ^ , 

for < t < 2tt. We compute that the maximum is achieved when t = tt/6, with 
the approximate value 416.269, and in particular 

IIL^II < 416.27. 

Substituting the estimate for the values of for i — 1,2,3, the estimate 

(5.11) for ||L(/)i||, and the estimate (5.12) for \\g\\ into the estimate for ||G|| ( |5.16 ), 
we arrive at our numerical estimate for \ \G\\ 

(5.18) ||G|| < 4.7772. 

We use this to estimate ||VG|| using the Poincare inequality 

||VG|| < yX^IIGII < W— 4.7772 < 69.32. 



Recalling the estimate for A (5.14), we have 

(5.19) |^|<JM|±MM) < t(1.243). 

V Vl-*^ll.9lP / 

Based on these estimates, we shall use the variational principle to estimate A2(i) 
from below. Since F — Acjji — (p + tG is orthogonal to the variational principle 
for A2 gives 

^ -J{F- A0i)Ao(F - A0i) _ -JiF- A0i)(A - tL){F - A0i) 

- l+i2||G|P • 

We compute the numerator to be 

\2{t)+t J FLF + J AM^ - tL){F - A(j)i) + J F{A ~ tL)A(t)i. 
The last two terms are 

A J <j)iAoiip + tG) + A J FAo0i. 

By definition of (f>i and integration by parts, these are 

= -AAi J M'P + tG) - J F(j}i. 

The first integral vanishes since (p and G are orthogonal to ipi. So, the numerator 
of ( [5:201 is 

- JiF~ Mi)^o{F - A(f,i) = A2(t) + t J FLF - AAi J Fc^i. 
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Thus, the variational principle for A2 implies 

^ \2{t)+tJ FLF-AXi J Fc^i 

i+i^iidp 

which gives the estimate for X2{t), 

\2{t) > X2il + f\\G\\^) ~ t [ FLF + AXi [ F(j)i. 



This implies 

A2(t) > A2 - t J FLF + AXi J F(j)i. 



Substituting (5.13), we have 

(5.21) X2{t) >X2-t J FLF + A^Xi > X2 - t J FLF. 

By definition of F = ip + A(f>i + tG and integration by parts, 

FLF = J LpLtp+A^ J (t)iL(j)i+t^ J GLG+2A J ipL(j)i+2t j GL^+2At / GLcj^i. 
This gives 

A2W >X2-tf (^(Lilt^o)^' 



-t J ip{L - Li\t=o)ip 
~tA^ I (PiL(j)i ~ I GLG - 2At / <fiL(l>i - 2t^ I GLip - 2At^ I GL(j>i. 



Incorporating estimate (5.1) for Ai(i), we have 

Ht) - Ai(i) >X2-Xi+t [ 0i(Li|t=o)0i - t I (^(ii|t=o)v - \Oi{t')\ 



-t j ip{L - Li\t=Q)(p 
-tA^ I (j)iL(j)i - I GLG - 2At / ^L^i - 2t^ I GLip - 2At^ I GL<j>i. 



Our calculations from the proof of Theorem 2 and the estimate (5.1) of |Oi(t^)| 
imply 

> ^ + (2.64)t - t2(i75.95) j ^{L~ ii|t=o)<^ 
2 / j, TA +3 I r<Tr< 0/1+ I ,„T A. 0+2 I r~<T,„ o /I+2 



~tA^ J 0iL</>i - r y GLG - 2At / ipL<j)i - 2t^ / GLip - 2Ar / GL(j>i. 

Recall the calculation 

4ta^9^ - 4t6^ag ^ -4^/3 , 4 ^ , ^ ^ , 

i - ii t=o = ^ ^ + — 7= 7= [aOxdy + bd 

' (\/3 + 2t)2 1(73 + 2^6)2 y 



y) 



Thus, we have / if{L — Li\t=Q)ip — 

Ata^ ,, ,|2 4t&2 ( 4^3 4 



ll'^o^ll + , , ^y \Vy\Vn , , - ^ U / <^x<^y+fo||<^,J| 



(V3 + 2t)2"^"" (V3 + 2i)2"^^" \(y3 + 2i6)2 ^3 
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We estimate using the Cauchy inequality, \\(px\\ and WfyW < ||V(/3|| with ||V(/j||' 
A2, 



w(L - Li\t-o)ifi < — 7= A2 + 2 — —-^ 



A2 



Since we assume t < 0.0004, we have 

J Lp{L - Li\t=o)ip < 0.58997. 
We estimate the remaining terms using the Cauchy inequahty, our estimates for 



G, IIVGII, and the general estimate (5.9 1 for L 



< \\Lcl)i\\ < 86.194, 



GLG 



< ^^ + ^^ ||VG|p< 22664. 



(V3 - 2ty 



ipL(t)i 



GLip 



GL<Pi 



< ||i0i|| < 86.194, 



< ||G||||L^|| < 2009, 



< ||G||||i0i|| <416, 



since \ \(t>i\\ = WfW = 1. Substituting the estimate for A gives our eventual estimate 
for the entire 0{t^) error term, and we have 

i{T{t)) > + (2.64)t - ^2(175.95) _ t(0.58997) 

-<3(1. 243)2(86. 194) - t3(22664) - 2*^(1. 243)(86. 194) - 2i2(2009) 
-2t3(1.243)(416). 



This becomes 



e(T(t)) > ^ + (2.05003)t 
9 



-i^ (175.95 + t * 133.174 + t * 22664 + 214.278 + 4018 + t * 1034.18) , 
Thus, we compute the largest t for which 

2.05003 >t{t* 23831.4 + 4408.23) . 
This is satisfied for any t < 0.0004. 



□ 



6. Proof of Theorem 2 

By our preceding results and continuity of the eigenvalues, we may now complete 
the proof of Theorem [3] by computing the first two eigenvalues of a large but finite 
number of triangles. 
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6.1. Continuity estimate. The following calculation is based on the linear defor- 
mation theory at the beginning of Section 4. We use T{x, y) to denote a triangle 
with vertices (0,0), (1,0), and {x,y), and we use Xi{x,y) to denote its i*'* Dirichlet 
eigenvalue, and £,{x,y) to denote its fundamental gap, 

^{x,y) = X2{x,y) - Xi{x,y). 

If a triangle T{x*,y*) satisfies 

{x* - xf + {y* - yf < t\ 



then by (4.1) 



2 M 

(6.1) i{x\y*) > ax,y) - (A2(x,y) + Ai(x,y)) . 

y 

Therefore, for each triangle T{x, y) at which we compute numerically 

^{x.y) > 



9 



we may use (6.1) to determine a neighborhood of triangles satisfying 

, 647r2 



9 ' 

without numerically computing the eigenvalues of the triangles in this neighbor- 
hood. Consequently, we have reduced the problem to numerically computing the 
fundamental gap of finitely many triangles and using the following algorithm. 

6.2. Algorithm. The main idea of the algorithm is to use the preceding calcula- 
tions to compute, to sufficient numerical accuracy, the first two eigenvalues of a 
finite grid of triangles and use this grid together with the continuity estimate to 
demonstrate that the gap of all triangles lying outside the cases covered by Proposi- 
tions 1 and 4 is strictly larger than that of the equilateral triangle. In particular, it 
follows from Propositions 1 and 4, the invariance of the gap function under scaling 
and symmetry that we need only compute for those triangles with vertices 

(0,0), (1,0) and (a;,y), 

such that the following inequalities hold. 

1.1 -|- < 1, by invariance of the gap function under scaling. 

1.2 0.5 < a; < 1, by symmetry. 

1.3 0.005 < y < 1, by Proposition 1. 



i.4 \J{x- 1/2)2 + (y - ^/3/2)2 > 0.0004, by Proposition 4. 

6.2.1. The steps of the algorithm. We begin with the triangle whose vertices are 
(0,0), (1,0) and (0.5,0.005); this is step 0. Next, in steps 1-2, we compute using 
( |6.1[) the radius t of the neighborhood around which the gap is strictly larger than 
647?/9. We then increase the x-coordinate in step 3, and check that the inequalities 
i.l-i.4 hold. If so, we repeat the calculations in steps 1-2 for the triangle whose 
third vertex is at the same height y but has been translated in the positive x- 
direction (to the right). We repeat steps 1-3 until the x coordinate is large enough 
so that one of the inequalities i.l-i.4 fails; then we proceed to step 4. In step 4, 
we return the ^-coordinate to 0.5 and increase the y-coordinate and check that the 
inequalities i.l-i.4 hold. We then continue repeating steps 1-4. 
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0. Initially, we define 

xo.j := 0.5 for all j, and yo := 0.005. 

1. At the {i,jy^ iteration of the algorithm, where the first iteration of the 
algorithm is (i, j) — (0, 0), for the triangle with vertices 

(0,0), (1,0) and {x,^j,yj) 

we compute 

1.1 Ai and A2, 

1.2 ^ = A2 - Ai and 

1.3 A = A2 + A1. 

2. We compute 



647r2 \ 



9 / 2.4^ 



Let riij be the smallest n e N such that the 10 " digit in the decimal 
expansion of t[ ^ is positive; let this digit be dij. Then, define 

t'i^j • — 10 '^di j. 

The numerical method must be accurate up to the 10^"^^ decimal place; it 
then follows from the continuity estimate (6.1) that for all triangles whose 
third vertex (x, y) lies strictly within a neighborhood of radius j about 
{xij,yj), the gap function is strictly larger than 
We define 

and verify the following inequalities. 
3-1 (x ^+i,,-)2 + (y^.)2 < 1. 

y(a;,+i,, -0.5)2+ (y,- > 0.0004. 

If these inequalities are satisfied, repeat steps 1-3. As soon as one of these 
inequalities is not satisfied, proceed to step 4 below. 
Let xo_j — 0.5, and for j > 1, define 



3.2 



and verify the following inequalities. 



4.1 (.Toj2 + (y^.)2 < 1. 



4.2 ^{x,^j - 0.5)2 - ^) > 0.0004. 

If one of these inequalities is not satisfied, then the algorithm is complete. If 
the inequalities are all satisfied, return to step 1 and repeat the algorithm. 

□ 



6.3. The numerical methods. The numerical computation of the eigenvalues 
were done by Timo Betcke using the Finite Element Method FreeFEM++ (8|. For 
efficiency, the calculations are made at each step but not stored, with the exception 
of toj which must be stored until it is replaced by ioj+i- To demonstrate the 
behavior of the gap function numerically, Timo plotted the logarithm of the gap 
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function in the figure below. The grid points are parametrized so that each grid 
point corresponds to a triangle with vertices (0,0), (1,0) and {x,y) where 

x^l-^ and y = ^^4- (2-r)2. 

Hence, the equilateral triangle corresponds to v — t — 1. 



11 




T 



Figure 2. Plot of the logarithm of the gap function on the moduli 
space of triangles 



6.4. Concluding remarks. Based on the numerics, we make the following con- 
jecture. 

Conjecture 1. The logarithm of the gap function on the moduli space of triangles 
is a strictly convex function. 



Recently, Laugesen and Siudeja 13 proved an interesting related result. 
Theorem 5 (Laugesen-Siudeja). For any triangle of diameter 1 with eigenvalues 

{Afc}r=i. 

n n 

(6.2) ^Afc>^Afc, V neN, 

k=l k=l 

where {^k}'kLi '^'"'^ eigenvalues of the equilateral triangle. 
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For n = 1, (6.2) is well known. The case n — 2 can be deduced from Tlieorem[3] 



as follows. By Theorem |3] and (6.2) with n — 1 



A2 - Ai > A2 - Ai =^ A2 + Ai > A2 - Ai + 2Ai = A2 + Ai. 

The existence and identity of a gap-minimizing simplex is a challenging open 
problem. Based on our results, we expect the following. 

Conjecture 2. Let DJln be the moduli space of all n-simplices with unit diameter. 
For n> 2, the regular simplex defined by points pQ,pi, . . . ,pn G M" such that 

\Pi - Pj\ = I for < i ^ j < n 

uniquely minimizes the gap function on 9Jt„ . 

There are several difficulties to be addressed. A subtle problem is the behavior 
of the gap of a family of collapsing simplices when several directions collapse simul- 
taneously. Is it possible that competing collapsing directions may result in a gap 
which stays bounded or converges to that of the interval as simplices collapse? Nu- 
merical calculations would provide insight into what one might expect; combining 
classical techniques with modern computation may produce interesting new results. 

We end this paper with a brief discussion of the similarities and differences be- 
tween the behavior of the gap function on convex domains and the gap function 



restricted to the moduli space of n-simplices. In the fundamental work of 16 



and subsequent papers 19 , 20 culminating in the proof of the fundamental gap 
conjecture the general method is to compare the eigenvalue estimate in higher 
dimensions to the eigenvalue estimate on a one dimensional manifold. The mini- 
mum gap for all convex domains can be asymptotically approached by thin tubular 
domains, and the minimum is achieved in dimension one. We pose the natural 
question: 

Is this minimum unique? 
More precisely, we make the following conjecture. 

Conjecture 3. Let D, C M" be a convex domain, and assume n > 1. Then 

In the case of triangular domains, the gap function is uniquely minimized by the 
equilateral triangle. It would be interesting to extend the beautiful works in the 
spirit of fl6' and [T to compare the eigenvalue estimate in higher dimensions to 
the eigenvalue estimate in dimensions greater than one. In particular, it would be 
interesting to compare the eigenvalue estimate to that on the equilateral triangle 
or other computable planar domains. 
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